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Abstract 

We determine the next-to-leading order renormalization group equations for the Two-Higgs- 
Doublet model with a softly broken Z 2 symmetry and CP conservation in the scalar potential. We 
use them to identify the parameter regions which are stable up to the Planck scale and find that in 
this case the quartic couplings of the Higgs potential cannot be larger than 1 in magnitude and that 
the absolute values of the iS-matrix eigenvalues cannot exceed 2.5 at the electroweak symmetry 
breaking scale. Interpreting the 125 GeV resonance as the light (7P-even Higgs eigenstate, we 
combine stability constraints, electroweak precision and flavour observables with the latest ATLAS 
and CMS data on Higgs signal strengths and heavy Higgs searches in global parameter fits to all 
four types of Z 2 symmetry. We quantify the maximal deviations from the alignment limit and find 
that in type H and Y the mass of the heavy (7P-even ( CP-odd) scalar cannot be smaller than 340 
GeV (360 GeV). Also, we pinpoint the physical parameter regions compatible with a stable scalar 
potential up to the Planck scale. Motivated by the question how natural a Higgs mass of 125 GeV 
can be in the context of a Two-Higgs-Doublet model, we also address the hierarchy problem and 
find that the Two-Higgs-Doublet model does not offer a perturbative solution to it beyond 5 TeV. 
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I. INTRODUCTION 


After the discovery of a scalar particle at the LHC [1, 2], one of the next qnestions is 
whether this is the one Higgs particle predicted by the Standard Model (SM) or whether 
there are more generations of SU(2) doublets, like it is the case for the fermions. In that 
sense, the simplest and most straightforward extension of the SM would be the addition of 
another Higgs doublet, the so-called Two-Higgs-Doublet model (2HDM) [3-5]. Furthermore, 
the measured mass of this new scalar [6] is a peculiar value for the SM: it tells us that the 
Higgs potential of this model cannot be stable up to very high energy scales [7, 8]. However, 
there is the possibility that the electroweak vacuum may just end up being metastable. So 
either one has to believe that we live in a metastable universe and then there is no need of 
new physics beyond the SM, or one has to introduce an additional mechanism to stabilize 
the Higgs potential. The latter could for instance be achieved by the heavier scalars of the 
2HDM. This model might be realized as an intermediate “effective” theory which describes 
physics at energy scales between the electroweak scale /igw of order 10^ GeV and some higher 
scale /ihigh- Beyond the latter, a more comprehensive model would be needed to describe 
“physics beyond the 2HDM”. An upper bound on //high is the Planck scale /ipi ~ 10^® GeV, 
at which gravitational effects become non-negligible in a quantum field theory framework. 
Large scale differences between //gw and //high bring along hierarchy problems like the fine- 
tuning of the 125 GeV Higgs mass, which could be resolved by mechanisms of the “complete” 
models, but are usually neglected in the effective models. Still one could ask to what extent 
the 2HDM could possibly mitigate the Higgs mass hierarchy problem and whether it might 
even be valid up to Planck scale without requiring any other New Physics. 

Therefore, we want to analyze the renormalization group evolution behaviour of the 
2HDM in this article, focussing on softly-broken Z 2 symmetric model realizations, which 
avoid flavour changing neutral currents at tree-level. Recently these models have attracted 
a lot of attention. A large number of papers [9-29] have analyzed current data for the 
125 GeV Higgs-like state within the context of 2HDM, and investigated the phenomenology 
of the other Higgs states present in the model. Given these results, the prospects for LHG 
upgrades and for other future colliders were examined in [30-37]. 

For renormalization group studies, especially the role of Higgs self-couplings is crucial 
and has been studied in the literature, in the SM (see for instance [7, 38]) as well as in the 
2HDM [5, 39-44], because these quartic couplings tend to destabilize the Higgs potential at 
some Phigh- Since a break-down of stability would mean that our theory would lose validity 
beyond a certain scale, we want to impose a stable Higgs potential beyond /Xew as a constraint 
on all couplings. Recently, the impact of stability up to the Planck scale on the parameters 
in the alignment limit of the 2HDM was discussed in [45]. If one wants to solve the Higgs 
mass fine-tuning problem, one has to guarantee the cancellation of quadratic divergencies 
of higher order Higgs mass correction terms. The corresponding conditions that need to be 
fulfilled are called “Veltman conditions” [46] in general, and in the context of the 2HDM also 
“Newton-Wu conditions” [47]. They have been analyzed at one-loop level [48-50] and even 
leading two-loop contributions have been taken into account in type H [51-53]. A recent 
idea was to only relax the cancellation of the generically large contributions of quadratic 
divergencies instead of imposing the strict cancellation using Veltman conditions [54]. 

In this article, we want to improve available results concerning two main aspects: We 
perform global parameter fits including the most up-to-date ATLAS and GMS results, rather 
than only using a handful of benchmark scenarios, which might not cover the whole spectrum 
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of interesting features. Secondly, we go beyond leading order precision by employing two- 
loop renormalization group equations (RGE) in order to analyze vacuum stability of the 
2HDM scalar potential. Moreover, we want to make use of the framework of next-to-leading 
order RGE to find out to what extent Veltman conditions can be fulfilled in the 2HDM. 

In the following, we first want to make the reader familiar with the model in Section II, 
and introduce in Section III its theoretical and experimental constraints and the numerical 
setup that we use. Then we are ready to compare leading and next-to-leading order renor¬ 
malization group equations for a benchmark scenario in Section IV A. We go on examining 
the quartic couplings varying the stability cut-off scale in global fits without experimental 
inputs in Section IV B. We also address the question of which upper limit to use for the 
unitarity condition from the perspective of renormalizability. Taking into account exper¬ 
imental data, we analyze the results of global fits to the physical parameters at pew and 
Ppi in Section IV C. The hierarchy problem is discussed in Section V, before we conclude 
in Section VI. Explicit expressions for the one-loop and two-loop RGE can be found in the 
Appendix. 


II. MODEL 


The Two-Higgs-Doublet model with a softly broken symmetry is characterized by the 
following scalar potential: 


+ A3($1$i)($^$2) + A4($I$2)(^2^i) + ^5 + (^2^l) 


( 1 ) 


where and $2 are the two Higgs doublets. In the following, we will use two sets of 
parameters: the eight potential parameters from Eq. (1), which we assume to be real, and 
the physical parameters consisting of the vacuum expectation value n, the GP-even Higgs 
masses rrih and uih, the GP-odd Higgs mass uia, the mass of the charged Higgs, mH+, the 
two diagonalization angles a and /3, and the soft Z 2 breaking parameter m^ 2 - The first two 
physical parameters can be treated as fixed by measurements, assuming that the 125 GeV 
scalar found at the LHG is the lighter GP-even Higgs. Instead of a and (5 we will use the 
combinations jS — a and tan P, since they can be directly related to physical observables. The 
measurements of the fight Higgs couplings to fermions and bosons are compatible with the 
SM, such that the 2HDM is pushed towards the so-called alignment limit [4, 31, 55, 56], in 
which P — a = TT/2. In this limit, it has recently been shown that CP violation in the scalar 
potential of Z 2 symmetric models with a soft breaking term is strongly suppressed [57], 
which qualifies our above assumption that the potential parameters are real. The masses 
of the heavy scalars could in general even be lighter than 125 GeV, and are not necessarily 
in the decoupling limit [4] (which itself is a limiting case of the alignment limit). In the 
following we will consider them to be in the range between 130 GeV and 10 TeV, that is 
heavier than the region where the 125 GeV scalar was found, yet still in the TeV range, 
which will be accessible by future colliders. 

Neglecting the first two generations of fermions, the Yukawa part of the 2HDM Lagrangian 
is 


Ty = - YtQpa2^*2tR - - l6,2(3L$2&fl - Yr^iL^^iTR - Yr^2LL^2TR + h.c. (2) 
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In the above Lagrangian, the top quark only couples to $2 by convention; its Yukawa 
coupling is related to the SM value hy Yt = / sin /?. Without breaking the Z 2 

symmetry in the Yukawa sector, there are only four possibilities to couple the Higgs fields to 
the bottom quark and tau lepton at the tree-level. They are called type I, type II, type X or 
“lepton specific” and type Y or “flipped”; in Table I we show the corresponding Higgs field 
assignments. Type H is of special interest, as it contains the Higgs part of supersymmetric 
models. As soon as we consider any one of the above types, only three Yukawa couplings 
remain as free parameters, and we can speak of 1^, W and Y- without any ambiguity. 


TABLE I. Yukawa assignments in the fonr possible Z 2 symmetric 2HDM types. 


Type I 

Type II 

Type X (“lepton specific”) 

Type Y (“flipped”) 

Y,1 = Yr,l = 0 
Y,2 = Y®^/sin/3 
Y2 = Y™/sin,5 

Yb,2 = Yr,2 = 0 

YM=Yfc®^/cos/3 

Yr,i^Y^^/cosl3 

Yb,l = Yr,2 = 0 

Y,2 = Y®^/sin/3 
Yr,i^Y^^/cos(3 

Yb,2 = Yr,l = 0 

Y,i = Y®^/cos/3 
Yr,2^Y^^/sm(3 


III. CONSTRAINTS AND SET-UP 

We will apply the following sets of constraints on the parameter space: On the theoretical 
side, the positivity of the Higgs potential [58] and the unitarity of the eigenvalues of the 
scattering matrix [59] are imposed at all scales and vacuum stability [60] at 
the electroweak scale. Moreover, we make sure that the quartic couplings A, (i = 1, 2,3,4,5) 
and the Yukawa couplings Yi {i = t, b, r) do not run into non-perturbative regions. On the 
experimental side, electroweak precision observables, the branching ratio Br(H Y^y), the 
mass difference in the Bg system and light and heavy Higgs searches constrain the 2HDM 
parameters at the electroweak scale. For a detailed description on the various constraints, 
we refer to [15] and [21] except for Br(i? —>■ Y^y) and the experimental Higgs data: In 
type H and Y we assume mH+ > 480 GeV in order to be consistent with the latest bound 
from Br(B —>■ Y^y) [61]. For the light Higgs signal strengths and heavy Higgs searches we 
use the most up-to-date ATLAS and CMS publications and pre-prints [62-79] , applying the 
narrow width approximation. We do not make use of (semi-)tauonic B decay observables, 
which would only be relevant in type H [80], because the existing tension between the 
measurements can only be explained in a 2HDM with explicitly broken Z 2 symmetry [81]. 
The SM parameters will be fixed to their best fit values [82]; for the SM Yukawa couplings 
in the MS renormalization scheme at the scale mz we take = 0.961, = 0.0172 and 

= 0.0102. While variations of the strong coupling aa{mz) within the 3cr allowed range 
have no effects on the outcome of our fits, varying the input for mt{mz) can have an impact 
on a specific parameter region like mentioned in [43]. However, we observe that these effects 
are imperceptible in the results of our global fits. 

The two-loop RGE have been obtained with the publicly available package PyR@TE 
[83]; we neglect all Yukawa couplings except for the top and bottom quarks and the r 
lepton. The observables have been calculated with the help of Zfitter [84-86], FeynArts 
[87], FormCalc [88], LoopTools [89], HDECAY [90-92], FeynRules [93] and MadGraph5 [94]. 
The frequentist fits are performed with the CKMfitter package [95]. For the fits involving 
experimental constraints we use of the naive definition of the p-value (Wilks’ theorem) [96]. 
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If not stated differently, exclusion limits are meant to be at the 2g level, which roughly 
corresponds to the 95% confidence level. 

Since we want to discuss various values for the scale /i in the following, we want to define 
our notation: The scale range for the running quantities lies between the electroweak scale 
/iew = at the lower end and the Planck scale /ipi = 10^® GeV at the upper end, as 
mentioned in the introduction. If - starting with a given set of parameters at /Xew and 
evolving to higher energy scales - one of the theoretical constraints is violated, we denote 
this breakdown of stability as /Xst- When discussing the hierarchy problem, it might be 
useful to introduce a cut-off scale \Xnat^ which a priori does not need to be the same as 
/ist- Furthermore, we want to introduce the parameter t(x) = ln(/i(js:)/ GeV) as the usual 
logarithmic scale. 


IV. RENORMALIZATION AT NEXT-TO-LEADING ORDER 

One can find a plethora of leading order RGE [5, 97-99] and next-to-leading order RGE 
[100, 101] for different realizations of a 2HDM in the literature; however, we failed to find a 
complete set for a 2HDM with soft breaking including the mass parameters, so we list the 
leading order (LO) and next-to-leading order (NLO) expressions in the Appendix. Before 
scanning over the whole parameter space with our fitting set-up, we want to explain some 
features of the 2HDM RGE looking at a representative example. 


A. A benchmark point 

In order to compare the LO and NLO RGE, we choose the scenario H-4 from [ 21 ] as 
benchmark scenario, because all quartic couplings are relatively large already at /Tew It is 
defined by mu = 600 GeV, vha = 658 GeV, = 591 GeV, P — a = 0.51371, tan /3 = 4.28, 
and mi 2 = 76900 GeV^, and compatible with all experimental measurements so far. The 
cut-off scale, where one of the quartic couplings becomes non-perturbative, is at 19.5 TeV at 
LO (dashed lines), and at 82 TeV at NLO (solid lines), see the top left panel of Fig. 1 . The 
Landau poles are at 54 TeV and 3.2 • 10® TeV, respectively; the former is shown as a vertical 
dotted fine in Figs. 1 and 2 . The fact that the higher order contributions “stabilize” the RG 
evolution, and thus increase both, cut-off and Landau pole scales, holds for all benchmark 
points from [ 21 ] and is a general feature in the 2 HDM: All dominant NLO contributions 
to the RGE of the A*, which are cubic in the quartic couplings, come with a negative 
coefficient and thus mitigate the positive LO contribution coming from quadratic A* terms 
(see Appendix). Beneath the total values of LO and NLO running we show the relative 
difference between LO and NLO RGE r* = |(AV — Af^°)/Af^°| with respect to the scale. 
For this benchmark point, the relative change of Ai and A 3 , is as large as 10% at around 
2 TeV and the difference increases even at a faster rate at higher scales. This is a first hint 
that the effect of the NLO contribution to the RGE in the 2 HDM is non-negligible. One 
can see that diverges around 35 TeV due to the fact that at this scale Ag^® turns to 0 . A 
better quantitative measure of the NLO vs. LO RGE is the relative distance 5^2 defined in 
[ 102 ]: Eor a dimensionless coupling L, we can define the relative distance between the LO 
and NLO curves (t) and L^^°(t) in the scale range from G to ^2 as 
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5 TeV 19.5 TeV 82 TeV 



FIG. 1. Leading order (dashed) and next-to-leading order (solid) RG running for the benchmark 
scenario H-4. On the left, we show on top the evolution of the quartic couplings. At LO, Ai hits 
the perturbativity limit dvr at 19.5 GeV; the Landau pole is at 54 GeV, indicated by the vertical 
dotted line. NLO RGE shift the perturbativity cut-off to 82 GeV. In the lower figure on the left, we 
show the relative error = |(A^^ — between LO and NLO expressions for the quartic 

couplings. The Yukawa couplings and the potential mass parameters are shown on the right. All 
types look the same except for the b and r Yukawa couplings. 


^12 — 


For H-4, 5^^ is 38%, if we integrate from mz to the LO cut-off at 19.5 TeV. To quantify the 
typical size of where Ac is the quartic coupling with the lowest perturbativity violating 
scale, we checked all benchmark points of [21], and found values between 17% and 45%, 
which indicates that in general the two-loop corrections are not negligible.^ 

In Fig. 1 we do not show the running of the gauge couplings gi, g 2 and 5 - 3 , since the 
two-loop corrections are too small to be visible. Also, the running of the Yukawa couplings 
is not significantly altered going from LO to NLO. However, due to the different assignment 
of the Higgs fields in the four types, we start with different values at the low scale (see 
Table I); that is why we denote the Yukawa couplings in the upper right panel of Fig. 1 as 
introduced in Eq. ( 2 ). (Note that only two of them are non-zero, depending on the type of 
Z 2 .) Among the mass parameters, m ^2 changes least if we run to higher scales, which we 
also observe as general feature of all types, mh and 771.22 have very different values at 
different scales, compare the lower right panel of Fig. 1. Neither of the mass couplings feeds 
back to the dimensionless couplings, as the RGE of the latter do not depend on m^ 2 ) ^11 


^ It is important to note that this definition of 6^2 is meaningful, if the denominator inside the square 
root does not become too small. 


tl _ 

tl 
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FIG. 2. Leading order (dashed) and next-to-leading order (solid) RG running of the physical angles 
and masses for the benchmark scenario H-4. The lines for m^+ and mn are almost on top of each 
other. At a scale of 1.4 TeV (2.8 TeV) at LO (NLO) v and ruh are 0 and /3 — o; is in the alignment 
limit of 7 r/ 2 . 


or m 22 - Furthermore, we have checked the mentioned benchmark scenarios for fixed point 
behavionr and do not find any below the perturbativity cut-off. 

If we switch to the physical parameter basis, we observe that also the RG running of 
the mixing angles can be sizable, see left side of Fig. 2. The scale at which (3 — a hits the 
alignment limit corresponds to vanishing v and rUh, which can be seen in the right panel 
of Fig. 2, where we show the running of all physical mass parameters. We find that the 
breakdown of the vacuum expectation value at some scale above fi^w is a general feature 
and occurs for all benchmark scenarios that we have analyzed; this is also observable in the 
benchmark points of [42]. 

After scrutinizing one benchmark point, we want to discnss more general features that 
can be fonnd in comprehensive fits. Especially the general dependence of /igt on the valne 
of \i{mz) will be an interesting question in the following. 


B. Fits without experimental data 

A parametrization independent way of setting upper limits to the quartic couplings of 
the Higgs potential is the requirement that the scattering matrix of $ 1 $^ $ 4 $^ processes 

is unitary. This corresponds to the condition that its absolute eigenvalues should be smaller 
than IGtt. The tree-level expressions [59, 103-106] are widely nsed theoretical constraints 
for the 2HDM; however, it seems that these bounds are very conservative. Studies involving 
higher order corrections have shown that the eigenvalues cannot be larger than 27r in the 
SM [38], and this bound has been adopted for the 2HDM of type II in [21]. Analyzing 
maximally allowed cut-off scales can shed light on how well this bound is motivated from 
the RGE perspective. 

In this section, we only want to impose the Higgs potential bounds, regardless of exper¬ 
imental constraints, in order to show the impact of the former on the 2HDM parameters. 
Since the assnmption of having a stable potential affects the potential parameters, we ex¬ 
press our results in terms of the five qnartic couplings and tan jS. (The latter modifies the 
Yukawa couplings as compared to their SM values, see Table I.) Due to the smallness of Yf, 
and Yr, their influence on the RGE is very weak and differences between the four Z 2 types 
are not visible in the A, planes. 
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FIG. 3. The blue (light) shaded regions show the dependence of Ust on the values for the quartic 
couplings Aj and tan (3 at mz if we choose an upper limit of 27r for the absolute S'-matrix eigenvalues. 
The red (dark) regions illustrate the allowed regions, if we take dvr instead and force at least one 
of the eigenvalues to be larger than 27r in magnitude. All types give the same dependence for the 
Aj. For tan/3, we show the possible regions in type I and X as shaded, and the areas below the 
dashed lines correspond to type II and Y. 


In Fig. 3 we show the dependence of the cut-off scale on the values of quartic couplings 
and tan/3 at the electroweak scale, for the two cases that either all eigenvalue moduli are 
smaller than 27r or that at least one of them is larger than 27r. Our fits show that forcing 
at least one eigenvalue of the ^-matrix to have an absolute value larger than 27r reduces 
the maximal cut-off scale pst to be at 5 • 10® GeV instead of the Planck scale; if we set at 
least one eigenvalue modulus larger than dyr, the maximal pst is at a few TeV. If we want 
to maintain a stable Higgs potential up to ppi, the largest eigenvalue can have a magnitude 
of at most 2.5 (~ O.Svr). Naturally, a larger upper bound on the eigenvalues allows for 
larger quartic couplings. But one can also see that cut-off scales larger than 10 TeV are 
only allowed for a very narrow range of tan/3 around 0.7 and - only in type II and Y - 
for an additional narrow range around 80. While the low tan /3 scenarios are known to be 
disfavoured for light 2HDM spectra by fiavour observables, we will see in the next section 
that also the large tan /3 regions are now excluded in type II. So we can conclude for all types 
but type Y that assuming Ust > 10 TeV and not too heavy new Higgs states all Y-matrix 
eigenvalues need to be smaller than 27r in magnitude. We will use the upper bound of 2?! in 
the following. Fig. 3 also shows the allowed \i{mz) and tan (3{mz) intervals for /igt at Planck 
scale. Roughly speaking, stability up to 10^® GeV requires |Aj(mz)| ^ 1 and tan/3(mz) > 1. 





































In this case, we also observe a lower limit on X 2 {'mz), which cannot be smaller than 0.15. 
In type II and Y, tan P{mz) is also limited from above and cannot be larger than 60. We 
list the precise ranges of the parameters in Table II. 

TABLE II. Allowed intervals for the quartic couplings and tan/3 at the electroweak scale, if we 
assume stability at pew (first line) and up to ppi (second line). 



\i{mz) 

X2{mz) 

Xz{mz) 

Mimz) 

X^{mz) 

tan/3(m^) 

Mst — Mew 

[0; 2.22] 

[0;2.20] 

[-1.8; 4.4] 

[-4.4; 3.8] 

[-2.1; 2.1] 

> 0.3 

Pst = Ppi 

[0;0.52] 

[0.15; 1.06] 

[-0.6; 0.8] 

[-0.9; 0.9] 

[-0.4; 0.4] 

> 1.0 in type I and X 
[1.0; 60] in type II and Y 


The bounds on \^{rnz) give us a handle on the question whether the Z 2 symmetry can 
be exact with stability up to the Planck scale: following [4], we find that the soft breaking 
parameter can be written as 


mi2 = 


tan P 


1 + tan^ (3 


[m\ + u^As) . 


( 3 ) 


Increasing /igt to higher scales not only gives a stronger lower bound on A 5 , but simultane¬ 
ously also excludes low rriA values, such that beyond pst ~ 10^° GeV a cancellation between 
the pseudoscalar and the A 5 contribution in (3) is no longer possible. Hence we confirm the 
LO result of [45] that a 2 HDM with pst > 10^° GeV has to be softly broken, which does not 
change significantly if we use NLO RGE. 

The inclusion of experimental bounds has only very little visible impact on the potential 
parameters, that is why in the following section we switch to the physical basis. 


C. Fits with experimental data 

In this section we want to show the impact of the experimental results discussed in Section 
III on the physical parameter space at the electroweak scale, once assuming a stable scalar 
potential at /igw and once for stability up to /ipi. We put special emphasis on the dependence 
of mass parameters on the relevant angles in order to investigate how large deviations from 
the alignment limit can still be. 

In Fig. 4, we show the tan /3-(/3 — a) plane for type I on the upper left, for type II on 
the upper right, for type X on the lower left and for type Y on the lower right. For a stable 
potential at the electroweak scale (orange) we show the la, 2a and 3a allowed regions (the 
2 cr region is shaded, the la and 3 cr contours are defined by the dash-dotted and dashed lines, 
respectively), and for a stable potential at Planck scale (purple shaded) we only present the 
2cr region. With stability at /iew, tan/3 is not constrained by any observable. For 2HDM 
masses below I TeV, however, we find a lower limit of 0.7 in all types (cf. [107]) as well as 
an upper limit of roughly 60 in type II. In contrast, /3 — a is constrained in all types to be 
fairly close to the alignment limit; the exact limits can be found in Table III. In type I, the 
deviations from /3 — a = 7 r /2 can be as large as O.Itt for a broad range of intermediate values 
of tan j3. Only a narrow band which is compatible with all constraints and at the same time 
allows for deviations from the alignment limit by more than O.OOtt survives the type X fits; 
within this band tan /3 is larger than 6 . In type II and Y, this band would in principle also 
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FIG. 4. tan/3-(/3 — a) plane in type I (top left), type II (top right), type X (bottom left) and 
type Y (bottom right) at mz with stability imposed at /iew in orange (light) and at /ipi in purple 
(dark). The dash-dotted, continuous and dashed lines border the Icr, 2a and Sa allowed regions, 
respectively; the 2a region - which roughly corresponds to the 95% C.L. area - is shaded. 


exist, but the new determination of the lower bound on mH+ from Br(S —?• Xg'j) excludes 
scenarios which feature 2HDM heavy scalar masses below 350 GeV and cut away the “lower 
branches” in the tan /3-{/3 — a) plane. This allows us to exclude a deviation by more than 
O.OStt from the alignment limit in those two types of Z 2 symmetry. We have seen in Section 
IV B that imposing stability up to /rpi constrains the quartic couplings; at this point, we 
want to shed light on the effect on the physical parameters. In Fig. 3, we already observed 
that tan /5 is constrained from below in type I and type X and additionally from above in 
type II and type Y. In type I we can also observe that for figt = /^pi, P — a has to be closer 
to 7r/2 for low and high values of tan/5 than in the case of figt = /igw In type X, the “lower 
branch” only occurs at 6.8 < tan /5 < 26 now, and also in type II and Y the allowed region 
is stronger constrained. Interestingly, the lower bound of 1 on tan/5 is not necessarily the 
same if we impose the alignment limit; in type II and Y we hnd tan /5 ^ 2 in this case. The 
reason why this value is smaller than the one found in [45] is that we use NLO RGE. At 
leading order, we conhrm their result that tan /5 < 3 is excluded in the alignment limit. 

In Fig. 5, we show the dependence of the charged Higgs mass bounds on tan /5; let us hrst 
discuss the case fist = In type I the strongest constraint for low tan/5 values comes from 
the mass difference in the Bg system. The other observables have no visible impact on this 
plane. The same holds for type X, except for mH+ < 300 GeV, where direct Higgs searches 
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TABLE III. Allowed intervals for /3 — a at the electroweak scale (and its sine and cosine) for all 
types of Z 2 symmetry, if we assnme stability at /Xew (first three lines) and up to /ipi (lines four to 
six). 




Type I 

Type II 

Type X 

Type Y 

Mst — Mew 

(3 — a 
cos(/3 — a) 
sin(/3 — a) 

[1.14; 1.91] 
[-0.33; 0.42] 
[0.908; 1] 

[1.49; 1.64] 
[-0.068; 0.081] 
[0.997; 1] 

[1.24; 1.70] 
[-0.13; 0.32] 
[0.946; 1] 

[1.50; 1.63] 
[-0.057; 0.076] 
[0.997; 1] 

Pst = Ppi 

(3 — a 
cos(/3 — a) 
sin(/3 — a) 

[1.21; 1.87] 
[-0.30; 0.36] 
[0.934; 1] 

[1.55; 1.62] 
[-0.044; 0.018] 
[0.999; 1] 

[1.29; 1.61] 
[-0.04; 0.27] 
[0.962; 1] 

[1.55; 1.61] 
[-0.040; 0.018] 
[0.999; 1] 


additionally cut away low tan/5 values. For type II and type Y, Br(i? Y^y) yields a lower 
limit of 480 GeV on m^f+; for large masses and low tan /5, the bound from the mass difference 
in the Bg system is stronger. In case of the type II we also find that a light charged Higgs 
is excluded for large tan P values; for instance if tan /? = 30, we obtain > 700 GeV. 
This is an effect only visible in a global fit: for large tan/3, heavy Higgs searches (mainly 
the tauonic decays) exclude light uih and m^. Electroweak precision data, however, are not 
compatible with too large mass splittings between the heavy neutral and the charged Higgs 
particles, so also the charged Higgs cannot be too light if tan /3 is large. This also qualifies 
that we did not use data from (semi-)tauonic B decays, which would give a weaker bound on 
the same corner of the type H plane. Type Y also features this constraint from the neutral 
Higgs searches, but it is much weaker and would only be visible for tan /3 > 100 because the 
T and b couplings to H and A cannot be enhanced simultaneously. Requiring stability up to 
/ipi gives almost the same regions as with stability at only that tan/3 gets constrained 
at the borders to stay within the limits from Table H. 

While the charged Higgs searches mainly depend on tan/3, neutral Higgs signals strongly 
depend on the deviation from the alignment limit, i.e. the actual value of /3 — a. Therefore, 
we show in Fig. 6 the allowed regions in the (/3 —a)-mj^ and (/3 —«)-m^ planes. For all types 
we observe that for neutral masses above 600 GeV the deviation of /3 — a from 7 r /2 can be 
O.OStt at most due to the stability bound. The larger deviations in type I and X correspond to 
neutral masses below 500 GeV, where the heavy Higgs searches become relevant constraints. 
As explained above, these regions are indirectly excluded by mH+ > 480 GeV in type H and 
Y and we obtain lower limits of 340 GeV and 360 GeV for uih and respectively. This 
lower bound on the pseudoscalar mass translates directly into a bound on the question 
whether the Z 2 can be exact, and combining Eq. (3) with the information from the allowed 
A 5 range in Fig. 3, we can conclude that even with a stability cut-off at the electroweak 
scale 771.^2 = 0 is very hard to achieve in type H and Y. If we additionally impose stability 
up to the Planck scale,we can see that sizeable deviations from the alignment limit are only 
possible for mn < 250 GeV and < 230 GeV in type I. Type X fits do not allow for /3 — a 
deviations larger than 0.027r for heavy neutral scalar masses above 150 GeV. In type H and 
Y, the lower bounds on the neutral masses increase to mn > 460 GeV and > 455 GeV, 
because in general, higher stability cut-off scales allow for less freedom in the mass splittings 
between mn, ttia and mH+ [9, 45]. In our fits we find an upper limit of 45 GeV on the 
absolute mass splittings for all Z 2 symmetry types. 
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FIG. 5. tan/3-mj:/+ plane in type I (top left), type II (top right), type X (bottom left) and type 
Y (bottom right) at mz with stability imposed at pew in orange (light) and at /ipi in purple 
(dark). The dash-dotted, continuous and dashed lines border the Icr, 2a and Sa allowed regions, 
respectively; the 2a region - which roughly corresponds to the 95% C.L. area - is shaded. 


V. THE HIERARCHY PROBLEM 

As we have already mentioned, there is a large scale difference between the Planck scale 
and the scale at which electroweak symmetry breaking occurs. This gap leads to the hierar¬ 
chy problem of the Higgs mass: if loop corrections can be of order of /rpi, why do they cancel 
each other almost perfectly, such that the Higgs mass is 17 orders of magnitude smaller? 
The cancellation of these mass corrections to retain a naturally light mu was first proposed 
by Veltman [46], therefore also referred to as “Veltman conditions”, and was first applied 
at leading order to the 2HDM by Newton and Wu [47]. Unlike in supersymmetry, in the 
2HDM there is no mechanism which naturally accounts for these cancellations. Still, there 
might be a hidden symmetry of which we are not aware, so nevertheless it is interesting to 
address this question. In the framework of the 2HDM, this hierarchy problem does not only 
affect nih but in principle also the other scalar masses, if they are not in the decoupling 
limit. However, since we do not know whether the heavier scalars are decoupled or not, 
we will only discuss the hierarchy problem of the already discovered 125 GeV scalar in the 
following. 

The largest one-loop contributions to the Higgs mass come from terms that are quadratic 
in /inat) if we assume that the 2HDM is valid up to a given scale /inat and use this scale 
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FIG. 6. (/3 — a)-mH plane (on the left) and (/3 — a)-mA plane (on the right) in type I, type II, 
type X and type Y (from top to bottom) at mz with stability imposed at /iew in orange (light) 
and at /ipi in purple (dark). The dash-dotted, continuous and dashed lines border the Icr, 2(j and 
3(7 allowed regions, respectively; the 2a region - which roughly corresponds to the 95% C.L. area 
- is shaded. 
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as a cut-off. Leading higher order contributions get an additional factor of [ln(/inat//^ew)]”, 
where n -|- 1 is the number of loops. If /inat is large enough, the logarithmic factor might 
compensate for the loop suppression, and the power series of the higher order corrections 
no longer converges. So requiring the cancellation of the first order Higgs mass correction - 
like often applied in the literature [40, 49, 50, 54] - is not sufficient if we do not know about 
the higher order terms. Only if we assume perturbativity of the power series, we can make 
a valid statement about whether the Higgs mass at the electroweak scale can be natural in 
the 2 HDM or at least whether the hierarchy problem can be mitigated. This assumption 
of perturbativity is analogous to the one applied above on the Yukawa and quartic Higgs 
couplings. 

All higher order leading logarithm mass corrections proportional to are given by 




oo / \ ^ 

(in^ ) 

n =0 ^ ' 


( 4 ) 


As described in [108] , especially for low cut-off scales the power series can be perturbative. 
However, we need to be careful to keep the leading logarithm sufficiently large with respect 
to the lower powers in the logarithm assuming that the leading logarithm gives the largest 
contribution. The leading coefficient function can be derived from the one-loop Higgs mass 
corrections and reads as 


3 3 3 3 9 

/o(A) fi)5'i) = “ 2 cos(2«)(Ai — A 2 ) + -Ai -|- -A 2 + 2 A 3 + A 4 + -^g\ + -^2 

- cos"(a) [eni + ] - siii"(a) [6^“, + . 

In order to easily obtain the leading logarithm contributions to higher orders, we use the 
recursive formula derived by Einhorn and Jones [109] , relating the coefficient functions fn+i 
to fn and the running of the couplings: 


/n+l 4 ^) Pi) . .. 'y ] o r Pi) 

n -|-1 oL 

This recursive relation is based on the following assumptions: the new theory has only 
one mass scale (m/j), and the logarithmic factor has to be large enough to suppress the terms 
with lower powers of logarithms. Two-loop effects on the Veltman condition have already 
been applied to the 2HDM of type H using this approach [51-53]; the authors found that 
the Higgs mass hierarchy problem can be ameliorated. 

An obvious choice of (i-aat as cut-off would be the breakdown of one of the stability constraints 
Usti so we will use it for the moment. In order to analyze whether we can make a statement 
on the Higgs mass naturalness in a 2 HDM which is based on a reliable perturbation series, 
we want to define kn as the ratio of the n-th correction term of Eq. (4) to the one of order 
n — 1 : 


k 


n 


9i) 

fn—1 9i) 


In^ 
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A^nat[T^V] 


FIG. 7. The allowed size of the one-loop coefficient of the Veltman condition series 

depends on the naturalness cut-off /inat- Without experimental constraints (light blue shaded), 
|/o(A^, can be as large as 75. With the inclusion of the measurements (orange shaded), it 

gets strongly constrained to be smaller than 6. More or less independently of taking into account 
experimental data, we obtain an upper bound on /inat of 5.3 TeV. The plane shows the type I fit, 
which agrees with the type X fit. In type II and Y, the maximal /inat is already at 3.7 TeV. 

Apart from the obvious logarithmic dependence on /inat, tho depend on the cut-off 
scale also indirectly: the latter determines which values for the A^ are allowed, see Fig. 3. 
Now we can re-write Eq. (4) as 


5ml 


^nat 

IGtt^ 


fo{Xi,Yi,gi) 


oo n 


1+X n 

n=l 1=1 


( 5 ) 


Only if we impose a small value of the leading order coefficient function and a sufficiently 
small number for ki, we can guarantee a perturbatively stable mitigation of the hierarchy 
problem of ruh-, also assuming that the for £ > 1 are not too large. Note that if we choose 
/o(Aj, hi, (yfj) to be exactly 0, ki diverges. 

If we constrain the first two factors ki and ^2 to be smaller than 1 in magnitude and that 
\5mW < we observe negative ki and ^2 in most cases, independently of the type of Z 2 
symmetry. This indicates that the series is alternating, which in turn means that - except 
for pathological scenarios - a suppression of the first two ki factors should be sufficient to 
make the series relatively robust with respect to pertubativity. Cutting the series in Eq. (5) 
after the second term (i.e. setting k^ = 0), we find that the maximal /Xnat is in the TeV 
range for all types, depending on the value we choose for fo{Xi,Yi, gi). The blue shaded 
region in Fig. 7 shows this dependence for type I, taking into account only Higgs potential 
constraints (as in Section IVB). There is a lot of freedom for fo{Xi,Yi, g^), which indicate 
large cancellations between the leading order contribution and higher order terms. This calls 
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into question our assumption that the series can be cut after the second term. The inclusion 
of the experimental results (cf. Section IV C, orange shaded in Fig. 7) limits the choice of 
/o(Aj, Yi,gi) to be of order 1, and thus presumably stabilizes the perturbative series. In both 
cases, however, the maximal /inat is at around 5.3 TeV for very small values of /o(Ai, V) 5'i)- 
While we obtain the same results for type X, the maximal /inat is even lower (3.7 TeV) in 
type II and Y which is a consequence of the much stronger constrained parameter space. 

Finally, one could also impose the perturbativity of the power series in Eq. (4) as con¬ 
straint and define /^nat as its breakdown scale if it is smaller than /igt- This, however, would 
not alter the maximal Huat, aor would it constrain the 2HDM parameters stronger than the 
conventional constraints. It would leave us with the question of what happens beyond the 
breakdown of perturbative naturalness already at a few TeV. 

To put it in a nutshell: softening the Higgs mass hierarchy problem is very difficult in 
the context of a perturbative 2HDM and can be achieved only for very low cut-off scales 
hnat- Nevertheless, this is an improvement of one order of magnitude as compared to the 
SM hierarchy problem and might hint at a more complete model beyond the 2HDM at TeV 
scales. 


VI. CONCLUSIONS 

We obtain the two-loop renormalization group equations for all four Z 2 symmetric types of 
the 2HDM using PyR@TE and show that in general, two-loop corrections to the leading order 
one-loop expressions should not be neglected. We then apply these equations to improve the 
predictions of renormalization group evolution of the coupling parameters, putting a special 
emphasis on the quartic couplings A*, which are usually prone to run into non-perturbative 
regions. The relative distance between the LO and NLO curves of the Aj can be as large as 
45%. The quadratic couplings mh and 771.22 from the Higgs potential can vary by an order 
of magnitude between the electroweak scale and the perturbativity cut-off, while 771^2 ^ 
general found to be rather stable under RG evolution. We do not observe any fixed point 
behaviour in the regions with a stable Higgs potential. 

Imposing positivity and perturbativity bounds at all scales and stability of the vacuum 
at the electroweak scale, the magnitudes of the A, at which give a stable Higgs potential 
up to the Planck scale are found to be typically below 1; we also find lower limits of 0.15 
for A 2 and of 1.0 for tan/5. We have checked that these results are the same in all types. 
In type H and type Y we additionally get an upper limit of 60 on tan/5 with stability up 
to /ipi. Moreover, we address the question of which upper limit for the eigenvalues of the 
tree-level —?• scattering matrix is appropriate and show that as soon as at least 

one of the eigenvalues exceeds 27r in magnitude, the maximal scale up to which the Higgs 
potential can be stable is 5 • 10® GeV. It even reduces to 10 TeV in all types if we assume 
1 < tan /5 < 60. Imposing stability up to fipy leads to an upper limit of 2.5 on the magnitude 
of the eigenvalues. 

Including latest results from the LHG as well as all other relevant experimental data, we 
show the result of our fits for all four types. We observe that deviations from the alignment 
limit strongly depend on the value of tan /5; the maximal deviation of /5 — a from 7r/2 is 0.43, 
0.08, 0.33 and 0.07 in type I, H, X and Y, respectively. (This corresponds to deviations of 
sin(/5 — a) from 1 of at most 0.092, 0.003, 0.054 and 0.003.) Taking the stability constraint 
up to /Upi, the bounds on /5 — a become even stronger and allow for deviations from 7r/2 of at 
most 0.36, 0.05, 0.28 and 0.04 in the respective types. The searches for heavy neutral Higgs 
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particles exclude a light charged Higgs boson for large tan /3 in type II and for very large 
tan (5 in type Y. In the — ct) planes it is visible that deviations from the alignment 

limit by more than O.Ohyr are possible only for tuh and rriA below 500 GeV in the types I 
and X. In type II and Y we obtain lower limits of 340 GeV and 360 GeV on tuh and rriA-, 
respectively This makes it very difficult to realize models with an unbroken Z 2 symmetry 
in these two types even if the stability cut-off is only at the electroweak scale. Demanding 
that the Higgs potential is stable up to the Planck scale, these mass limits are even stronger. 

We finally discuss whether a reliable statement on the seemingly fine-tuned Higgs mass 
can be made in the context of a 2HDM and whether its hierarchy problem can be solved at 
least partially. Restricting higher order corrections to the perturbative regime, we observe 
a maximal naturalness cut-off at 5.3 TeV. Our conclusion is that within a perturbative 
framework a natural cancellation of quadratic divergencies cannot be implemented into a 
Two-Higgs-Doublet model beyond (P(TeV) scales. 
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APPENDIX 


Here we list the renormalization group equations for the 2HDM with soft Z 2 breaking 
which we obtained with the PyR@TE code [83] . 

Eor any coupling L the complete (3 functions at NLO can be split into leading and next-to- 
leading order contributions and further divided into bosonic and fermionic parts as follows: 


Pl — 


dL 

dt 


= 


l^(N)LO ^ ^(N)LO,b 


Except for the Yukawa RGE the bosonic part does not involve fermionic couplings and is 
type independent, while the fermionic part in general depends on the type of Z 2 symmetry. 
If the expressions for the latter differ for the different types, we will replace the index / by 
the type label I, II, X or Y, respectively. 
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The RGE of the gauge couplings only depend on themselves and on the Yukawa couplings: 




= - Sjf 

(1671^/3^^^^’’’ = {2gl + Sgl + 12 ( 72 ) gl 

= - ( 5 V- + 


1633^/?“ = - 7j| 

(1633^)^"“’^ = - (2n" + 21?) 93 “ 


As already mentioned, the mass parameters from the Higgs potential do not influence 
the running of the dimensionless couplings. Their running, however, is not negligible and is 
given by the following expressions: 


~ ~ 2^2 (^"^3 + 2A4) m22 

167r2;d^?’^ =0 

167r2/3^^’" = (6^2 + 2Y,2) ^ 2 ^ 

IQTT^P^f =2T>?1 
167r2/5^«;^ =6n2m2^ 
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Finally, the quartic couplings from the Higgs potential: 
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=(167r2)2/3™’" - \glY,%^ - %AglY^Y^ + 36^^^/ + 36^^^/ 

+ Y^Y^\m\^ + 24 A 4 ) - ( 3^2 + Y^) (Ai - A2)(12A3 + 4 A 4 ) 
=(167r2)2/5™ + y,2(Ai - A2)(12A3 + 4 A 4 ) 

IGtt^/?™’^ =(167r2)24™’^^ + Y^{\^ - A2)(12A3 + 4 A 4 ) 


l 67 r^/ 3 ^^^’^ —“ig^g^ — i^gl + ^gVj + 2A1A4 + 2A2A4 + 8A3A4 + 4A4 + 8A5 
levr^/^^f= (QY^ + QY^ + 2Y^) A 4 
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167r^/5^f =12Y^Y^^ + + 6Y^^ + 2Y^) A4 

1677^ =1Qt7‘^P^°'^ 

167r2/5^f=167r2/5^f 

f-\r 2\2nNLO,b '^^4 2 t a 2 4 

(167r ) ’ = - —g^g^ - Ug^g^ 

51 \ 231 

5Ai + 5A2 + 2A3 + -^fl'2'^4 

+ (^AXiX^ + 4 A 2 A 4 + 4 A 3 A 4 + 8 A 4 + IGA^^ + ^2 ( 36 A 3 A 4 + I 8 A 4 + 54A5^ 

— (TA^ + 7A2) A4 — (Ai + A2) (4OA3A4 + 2OA4 + 24A5) 

- 28A^A4 - 28A3A^ - 48A3A^ - 26A4A^ 

(16x2)"/?™’' (gj.2 g. 21 J .2 g. iiy^2) g. ^gj.2 g. ®y^2 g. 

+ A (f n" + f 1)2 + f 1 ?) A4 + 409" (If + if) A4 
-f (1)2-1)2)2 A 4 -?ifA 4 

— (121^^ + 121^^ + 41 ^^^ (A2A4 + 2A3A4 + A4 + 2A5^ 

(167r2)2/5f/^’^^ =(167r2)2/5f/^’^ + \glY,%^ + UglY.^Y^ - 2AY,%^ - 24^^/ 

- n'^i'(24A3 + 6OA4) - (12^2 + 4Y^){X, - A2)A4 

=(167r2)2/5™’^^ + 4y,2(Ai - A2)A4 


157 4, 2 2f 

+ + 5'lfi'2 ( 


IGtt'^P^^'^ — (— 3gi — 95^2 Y 2X1 Y 2X2 + 8A3 -|- I2A4) A5 
167r2/3^f’7 = {6Y^^ + 6 F/ + 2Y^) A 5 

+ ^gfgl - ^^2 - gf ( 2 Ai + 2 A 2 - 16 A 3 - 24 A 4 ) 

Y 5^2 (36A3 -|- 72A4) — 7A^ — 7A2 — (Ai -|- A2) (40A3 -|- 44A4) 




- 28A2 - 76 A 3 A 4 - 32A2 + 6A2 1 A, 


2 I 25 2 35 2 25 2 



+ fi '2 


45 2 35 2 15-4^2 

jlf + ^1)2 + -n" 


+ 932 (40if + 40yi") - 5 (if - y;^)2 - lif 


~ (l 23 ft^ + 12A^^ + YY^^ (A2 -|- 2A3 -|- 3A4) J A5 
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- {36Y,%^ + 12^2 (Ai - A 2 ) + 4F,'(Ai - A 2 )) A 5 

167r2/5™’^ =(167r2)2/3f^^04 _ 

167r^/5™’'^ =(167r^)^/5j^^’^^ + 4y,^(Ai - A2)A5 


Strictly speaking (and according to the introduced notation), there are only fermonic 
contributions to the (3 functions of the Yukawa couplings. However, we will denote the 
bosonic loop contributions which are the same in all types as Note that one 

bosonic contribution of and depends on the type; thus we add the term to 

the corresponding Py,^^'^. 


16?:^Py’^ = - {^gl + ^gl + Yt 

167r^/5^f+ ^Y/ + Y,^^ Y, 

IGrr^Py^’^^ = - {Y^ + Y^) Yt 

IQir^Py'^ = IQir^Py'^ - Y^Yt 

- Y^Yt 


(167r2)2/3^^o,fe ^ 


/1267 4 322, 22 ^^4 


- 4^i^ 2 + yPiP3 - ^^2 + 9^2^3 - 108^: 

3 


+ 2-^2 + -^3 + ^ 3^4 + A4 + -A5 ) Yt 


= )?(+ j^'r \+si [^ 



2 , 99^,J , 225^,J , 15^,2 


+ al (iy^ + 36Vf) 


2\2 nNLOJI _ 




IGtt /5 


2nNLO,X _ 


-p'i - jifVf + - I 2 y,- - - 5v;* - 6 V-A 21 r, 

2 , 41 , . 131 A 2 /^33,,2 225, A 2 /^16..2 0...2 

»■ ( -TTY + l6^‘ ) + (l6'‘ + 16^‘ ) + 

-pi* - p?y^ - p'^y^ - i2y,‘ - anUa+2v;U4 - 1 y, 

2\2nNLO,I /25___2 , 13^2 , ^^^2 ^ ^^2 ^ ^,^2 \ ^,^2^ 



IQTT'^Py^ 


Yt 


2 pNLO,Y 




2\2pNLO,II 


= i^Q^ypYt 


+ 


25 9 15 9 39 99 99x9 

—g^ + —gl + -Yh - -y? - -y^ yXf 

8 ^i 8^2^4& 4t 4 -'-^ 


l % T :‘^ ppy ^ = - (^gl + + 8^2^ Yb 
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+ \ y ? + if) n 

= IGt :^/ 3 y ^’^ - Y ^ Y ^ 
lQTT ‘^ f 3 ^^'^ = - Y^Yk 

leTT^/^ff= i 67 r 2 / 5 ^f- (y ,2 + Y ^) n 


/ 1 /- 2\2nNLO,b I 4 

( 16 ir ) • = (^-^ 9 , - jSife + J 


2\2oNLO,I 




2\2 pNLOJI 




(,2 [79 2 m ^ \ ^ ^ ' 

[ 9 i (^ign + 48 ^ + ^ 16 ^* + ° 

+ gl ( 36^2 + 4 y, 2 ) _ i 2 Y ,^ - ^ Y ,%^ - ^ Y,X 


+A3 + A3A4 + A4 + 2-^5^ ^ 


t: t: 

-iv;‘ + 5 f ,= v ;^ - ?y;.- - 6 Y ,^\, + 5 a =) n 


/ 2 / 79^^2 53 ^^2 25 ^^ 2 ^ 2 ^ 225^^2 33^^2 

(»■ W ^ - uY + + 16 ^- + ) 

+9I (sen" + f >;") - i 2 n‘ - jn'n" - jn'n" 

-^if - jif - 6 if Ai - 2 Yf \, + 2 r ,% + ? a ;) n 

( 16 ^")"/)“°s - (|j; + - ?F/ + 5 f ," - ?F,") F,"n. 


^ 225^^2 33^^2 15,^2' 

w + Te^-+ T^-, 


1 ^ 2 nNLO,X 
lOTT Py ^ ’ 




- 2\2nNLOJI 

ITT /■),. 


( 167 r^)V’ 


^ 9 5 2 _9 2\ 2 ^ 

U®' 8®" 4 ‘+4 ' 4 4 

/25 , 15 , 9,,2 , 3 2 9y2W2y 


2 nLO,b / 15 2 , 9 2 1 -IF 

IGtt = - I —+ -c /2 1 

167r2/5^f= 1^3^^," + 3y," + ^Y^^^ Y^ 

IGyr^^yf - SY^X 
IQtt ^ Py ?’^ = IQtt ^ Py ^’^ - pY ^ + 3 ^^,^) Yr 
= Wtt '^ Py ^’^ - YY^Yr 

f-xr 2\2QNLO,b / 161 4 922 21 4 .2 
(16vr ) Py^ ’ = ( —Pi + -P 1 P 2 - -j92 + ^3 
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+ si (20V? + 20V?) - ^y? + 2v?v? - ^v?v? 


- jv ; - jV ? y ? - 3 y ? - 6 y?Aj + ^ Ai ) vy 


2\2nNLO,II 


{■\r 

( 167 r ) 


25, 


179, 


24 


16 


= Ur TTtU + -TTrU ]+9n + -TTrU + ‘^^gtYk 


45, 


165, 


,2-y^2 


167r^/3' 


Vr 


= (167r2)2/5- 


Yr 


16 

~Y* - ?y?y? - i^y?y? - ay? - 6y?A: + ?a?) v; 

25 , 45 n n 27 , 9 , 27 n\ , 

+ y^2^ + 20,1 - -n^ - -U - y U 




, 25 2 45 2 ^. 2 27,^2 15..2 27 ,,., , 

- ( + yP2 + 20,2 - -y,2 + _y^2 _ _Y^. j Y^^Yr 


= ( 167 r 2 ) 2 / 3,^^04 / , -.2 , on.2 - X .2 , -^2 - ^a2 \ ^.2 
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